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Abstract. We show that a certain class of vortex blob approximations for ideal hydro- 
dynamics in two dimensions can be rigorously understood as solutions to the equations of 
Cu ' second-grade non-Newtonian fluids with zero viscosity, and initial data in the space of Radon 

^ , measures M.{M?). The solutions of this regularized PDE, also known as the averaged Euler 

or Euler-a equations, are geodesies on the volume preserving diffeomorphism group with 
respect to a new weak right invariant metric. We prove global existence of unique weak 
^^ ' solutions (geodesies) for initial vorticity in M.{M?) such as point-vortex data, and show that 

H , the associated coadjoint orbit is preserved by the flow. Moreover, solutions of this particular 

vortex blob method converge to solutions of the Euler equations with bounded initial vortic- 
ity, provided that the initial data is approximated weakly in measure, and the total variation 
of the approximation also converges. In particular, this includes grid-based approximation 
schemes of the type that are usually used for vortex methods. 
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1. Introduction 

The starting point of our investigation is the somewhat surprising fact that 
the equations of motion for an inviscid non-Newtonian fluid of second grade, 
and Chorin's vortex blob algorithm with a particular choice of cut-off or blob 
function are, at least formally, equivalent. 

The velocity field u = u{x, t) of a second grade fluid, under the assumptions 
of observer objectivity and material frame-indifference, satisfies the unique 
equation 

(1 - a'^A)dtU + u ■ V(l - o?l\)u - o? (Vm)* ■ Am = - gradp , (1.1a) 

divM = 0, (1.1b) 

m(0) = mo, (1.1c) 

where p = p(x, t) is the pressure function which is determined (modulo con- 



stants) by the velocity field. See ]TT| and references therein for a discussion of 
the constitutive theory of second grade fluids, and 0, ^ for well-posedness of 
the viscous second-grade fluid equations. In this context, the constant a > 
is a material parameter which represents the elastic response of the fluid. 
In two dimensions, taking the curl of equation ( |l.la ) and setting g = (1 — 



a^A) curl2D u yields the vorticity form 

dtq + u ■ grad g = , (1.2a) 

u = K''*q, (1.2b) 

g(0) = go, (1.2c) 

where q = q{x, t) is called the potential vorticity^ and K"" is the integral kernel 
of the inverse of (1 — Q;^A)curl2D, defined so that the divergence condition 



l.lbl) is satisfied. 



When a is interpreted as a length scale, (|1 . 1|) or (|1.2|) are known as the aver 



aged Euler or Euler-a equations [0 which model the large scale flow (spatial 
scales larger than a) of an ideal incompressible fluid. Their analysis and rich 
geometry has recently received much attention |^, |2^, |2^. In particular, so- 



lutions of (|1 . 1|) on an n-dimensional Riemannian manifolds (M, g) arise as 
geodesic flow on the group of ilf ^-class volume preserving diffeomorphisms T>1^ 
provided s > {n/2) + 1 with respect to a new weak right invariant metric, 
given at the identity element e G Vf^ by 

{u,v)e = {u,v)j2 + 2a^(Def u, Def f)^2 (1.3) 
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where Def m = (Vu + Vm*)/2. Thus, following the program of Arnold and 
Ebin-Marsden [|13| , local- in-time well-posedness of classical solutions is a direct 
consequence of the existence of C°° geodesies of ( ■ , ■ ) on P* . 

The vortex blob method was introduced by Chorin as a regularization of 
the point vortex algorithm for ideal hydrodynamics, and can be understood as 
follows. Consider the vorticity form of the Euler equations on R^, 

dfUJ + u ■ gradu = , (1.4a) 

u = K*uj, (1.4b) 

uj{0) = uJo. (1.4c) 

Here K{x, y) = l/(27r) V"*" log \x — y\ and u = uj{x, t) is the physical vorticity 
of the flow. When the velocity field is sufficiently regular — u is at least contin- 
uous in t and quasi-Lipschitz in x, uniformly over finite intervals of time — we 
may define the Lagrangian fiow map rjt = r]{- ,t) hj 

dtr]{x,t) = u{r]{x,t),t), (1.5) 

or equivalently by 

dtTjt = Ut O Tjt . (1.6) 

For each t, the map f]{- ,t) is in Q, the group of all homeomorphisms of M^ 
which preserve the Lebesgue measure. The pointwise conservation of vorticity 
under the Euler fiow is thus expressed by Utorjt = Uq; combining (|1.5|), ( |1.4b[ ), 
and the initial condition ri{- ,0) = e, we obtain the ODE 

dtvix, t)= f K{r]ix, t),ri{y, t)) uj^{y) dy . (1.7) 

Letting 6 denote the Dirac measure and substituting the point vortex ansatz 

N 

ujix,t) = Y,^i^ix-Xiit)), (1-8) 

n=l 

into ( |1.7|) , we obtain a finite dimensional system of ordinary differential equa- 
tions for the vortex centers Xi, . . . ,xn- However, the induced velocity field has 
l/|a;|-type singularities at the vortex centers. Hence, the point vortex system 
is neither numerically well-behaved (the exact solution of the point vortex sys- 
tem may even collapse in a finite time for small sets of initial data ||2^), nor 
does it approximate physically relevant velocity fields very well. 

The idea of the vortex blob method is to smooth the Dirac measure by a 
cut-off or blob function x that decays at infinity and whose mass is mostly 
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supported in a disc of diameter a. This leads to the following equation for the 
Lagrangian flow: 

diTr{x,t)= [ K^{r{x,t),r^^{y,t))u,{y)dy, (1.9a) 



where 

K" = V^G" , (1.9b) 

-AG"(x,y) =x"(|x-i/|) = Ixf^^l • (1.9c) 

a^ \ a / 

Many researchers have investigated the convergence properties of this scheme 
0, H, ^ ^ 1^. In particular, for certain smooth cut-off functions, such as 



Bessel functions, the order of accuracy with respect to the regularization pa- 
rameter a depends only on the smoothness of the Euler flow ( "inflnite order 
accuracy" ) . 

It is now easy to see that the equation of a second grade fluid ( |1.2| ) and the 
vortex blob method coincide when x{^) = Kq{x), where Kq is the modifled 
Bessel function of the second kind which is the Green's kernel for the operator 
(1 — A). Thus far, this relationship has only been formally established, as it 
remains to be proven that the point-vortex ansatz ( |1.8| ) makes sense as data for 
the PDE ( |1.2| ); moreover, it is not a priori clear if solutions to the vortex blob 
method converge to true Euler solutions (in the sense of PDE). Our results are 
the following. 

We show that the Lagrangian flow formulation of the blob method ( |1.9| ) with 
Kq cut-off function is well-posed for initial potential vorticities go in A^(M^), 
the space of Radon measures on M^. In particular, this includes point-vortex 
initial data. Such a result does not hold for the Euler equations, where the 
flow map of the point vortex system ( [1.7|) is not known to be well-deflned. 

This result allows us to rigorously classify the co-adjoint orbits characterized 
by point- vortex initial data. Let us explain what we mean by this. The conflg- 
uration space for ideal incompressible hydrodynamics is the volume-preserving 
diffeomorphism group, and for s > {n/2) + 1, "D^ is a C°° Hilbert manifold 
and a smooth topological group. While P^ is not a Lie group (left composi- 
tion and inversion are only C^ and the group exponential map does not cover 
a neighborhood of the identity), it behaves similar to a Lie group, because 
of the smooth properties of the Riemannian exponential map (see |1^ and 



231 p4|)- '^^^ Eulerian phase space for the fluid motion is the single flber 
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TeD^ consisting of if ''-class divergence free vector fields on the fluid container, 
and this vector space can be formally thought of as the "Lie algebra" of P^. 
The cotangent space at the identity is given by H'^ (A^) / dH'^ (A^) , the if ''-class 
differential 1-forms modulo exact 0-forms. Using the fact that the exterior 
derivative d: T*Vf^ —>■ ii^^^(A^) is an isomorphism, and the fact that we may 
identify H'-^^A"^) with H^-^A'^), the role of the dual of the "Lie algebra" for 
2D hydrodynamics is played by the ii''" ^-class vorticity functions. The rep- 
resentation of P^ on this "Lie algebra" is provided by the co-Adjoint action, 
so that for rj G P^ and u G if^~^(A°), Ad*(c<j) = 0007], and the invariance of 
the CO- Adjoint orbit is merely the pointwise conservation of vorticity which is 
fundamental to 2D hydrodynamics. If one temporarily ignores the topology 
and works formally, then it is possible to classify certain interesting and im- 
portant co-Adjoint orbits. Specifically, it is a result of Marsden and Weinstein 
p2| that point-vortex initial data ( [1.8|) define the co-Adjoint orbit on which 
point- vortex dynamics evolve. This is clearly a formal result as Dirac measures 
are not elements of H'^"^ for s > 2; consequently, the problem is to supply a 
candidate topology for the "Lie algebra" which is general enough to contain the 
Dirac measures, and weaken the regularity of the configuration space so that 
its "representation" is well-defined. In doing so, one can establish a rigorous 
classification of the orbit. By using Q for the configuration space and A^(M^) 
for the "Lie algebra," and by defining a new notion of weak co-adjoint action 
which coincides with the notion of a weak solution, we are able to establish the 
orbit classification for point-vortex initial data, and prove that our particular 
vortex blob method leaves such weak co-adjoint orbits invariant. 

Finally, we consider the matter of greatest practical importance: the conver- 
gence of solutions of the vortex blob method to solutions of the Euler equations 
as the blob diameter a ^ 0. We prove this convergence result under the rather 
mild assumption that the initial Euler vorticity field ujq is continuous with 
compact support and is approximated on its support by a sequence of weakly 
converging measures in A^(R^) that have uniformly bounded total variation. 
(The restriction to compact support will be replaced by weaker assumptions 
on the decay at infinity.) 

The precise statements of our results are as follows. 
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Theorem 1. For initial data go G A^(M^), there exists a unique global weak 
solution to (11.21) with 



ri" eC\R;g), M"GC°(M;C;Ji,(M^M2)), and g G C°(M; A^(M2)) , 

(1.10) 

where the subscript div denotes divergence- free. As a consequence, the co- 
Adjoint action Ad*(g) and the weak co-adjoint action w-ad*(g) are conserved. 

Remark 1. The solution that we construct may not necessarily have finite 
energy, i.e., the velocity field m" may not be in L^. None of our results, however, 
relies on energy type estimates. Furthermore, as is the case for the Euler 
equations, the initial potential vorticity can be decomposed into a radially 
symmetric and a mean-zero part, with a corresponding velocity field u in the 
affine space Mgtationary + i^^(M^,]R^). For details see DiPerna and Majda |T^. 

Remark 2. An immediate consequence of the uniqueness of the solution and 
the time-reversibility of the equation is that the vortex blob system cannot 
collapse in finite time, i.e., two or more vortex centers cannot merge into one 
in finite time. For non-regularized Euler vortex dynamics, on the other hand, 
it is known that vortex collaps occurs on small sets of initial configurations 

m. 



Remark 3. The kernel K"' which corresponds to a second grade fiuid is the 
least regular kernel (modulo possible sub-logarithmic corrections) for which 
uniqueness of point vortex solutions can be shown. An equivalent uniqueness 
result based on Sobolev space methods, and for bounded domains is given in 



15 



Theorem 2. Let r] be the flow map of the Euler equation ( |1.7| ) with initial 
vorticity luq G L^(M^)nL°°(M^). Suppose thatujQ is approximated by a sequence 
of measures q^ in M.{E?) such that q^ -^ ujq weakly in M.{E?) and ||gollA^ ~^ 
llciJoll^^i. Then for every T > 0, there exists a sequence {an} converging to zero 
as n —y oo such that when rj"'" denotes the flow map of the vortex method with 
a = an and initial data q^, 

lim sup sup [//"""(x, t) — r7(x, t)| = . (l-H) 

"^°^ie[o,T]2^eiR2 
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Remark 4. The idea of analyzing the vortex method as a PDE posed on some 
space of distributions was aheady used by Marchioro-Pulvirenti |]19[ and Cottet 



0. Cottet 's result requires stronger assumptions on the cut-off function and 
hence smoothing kernel, stronger regularity assumptions on the underlying 
Euler flow, and his approximation of the Euler initial vorticity field required a 
uniform grid. The trade-off, however, is that these more stringent constraints 
give an improved (algebraic) convergence in a. There is, in general, a trade- 
off between the order of convergence on the one hand, and the assumptions 
placed on K"", cjq, and the approximation at time t = on the other. The result 
which, to our knowledge, comes closest to Theorem ^ is given in Marchioro 



and Pulvirenti |19|. The authors, however, assume that K°' is Lipschitz, which 



again excludes kernels corresponding to the equations of second grade fluids. 

Remark 5. Stronger results can be proved for kernels K"" with a higher degree 
of smoothing. For example, by replacing 1 — a^A with (1 — a^A)'', one obtains 
a hierarchy of regularizations of the Euler equations which coincide with geo- 
desic flow on the volume-preserving diffeomorphism group with respect to the 
H^ metric. Other choices of K"" may introduce non-local pseudo-differential 
operators into equation (|1.2|), but the analysis can still proceed as before. 



Remark 6. In three dimensions, the formal connection between second grade 
fluids and particular vortex filament methods still holds and is the subject 
of a forthcoming article. In this setting, one looks at the set of vorticity 
distributions of the following form. Let 7 be a curve in M? extending to infinity 
in both directions, and let 5^ be the Dirac distribution given by integration 
along 7 with respect to arc length. Let uj^ be the 2-form along 7 defined by 
ixdz A dz, where T is the unit tangent vector to 7. Then if F is any constant, 
Voj^d^ is the vorticity corresponding to 7 with strength F. See WI. 



Remark 7. As we described above, for s > {n/2) + 1, local well-posedness 
follows form the existence of unique C°° geodesies 'r]{t) on V'l with respect to 
the right invariant metric (■, ■) defined in ( [I.3[ ), with initial conditions 77(0) = 
e and 77(0) = mq. In working with the geodesic flow f]{t), one obtains C°° 
evolution in the tangent bundle TP* and C°° dependence on initial data, 
while the projected evolution curve u{t) = fi{t)or]{t)~^ in the single fiber of the 
tangent bundle T^/D^^ - which plays the role of the Eulerian phase space - has 
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only C° smoothness, and C^ dependence on the initial velocity field. In the case 
that the manifold has a smooth boundary, there are three new subgroups of V^^ 
which are in one-to-one correspondence with the classical Dirichlet, Neumann, 
and mixed elliptic boundary value problems in the sense that elements of the 
"Lie algebras" of these three subgroups satisfy those boundary conditions. 
Hence, geodesic fiow of (■, ■) on these three subgroups gives the solutions of 
( |1 . 1| ) with no-slip, free-slip, and mixed boundary conditions. 

2. Kernel estimates 

The crucial ingredients for the proof of our theorems are quasi-Lipschitz 
estimates on the Euler kernel K and the regularized kernel K"". For x G M^ 
we define the function 

^^^) ^ JN (1- In N) for|x|<l, ^2.1) 

1 for X > 1 . 



Lemma 3. For uj G L^{m?) n L°°(R2), 

\K{x,y) - K{x\ y)\ \uj{y)\dy < c^{x - x') {\\uj\\^, + ||^||^oo) . (2.2) 

The proof is standard and can be found, for example, in McGrath ||18[|. 



Somewhat less standard is the following estimate, still for the Euler kernel, 
which is similar to estimates in Benedetto et al. [H. 

Lemma 4. Let uj G L^(M^) fl L°°(M^) and let (p he an area preserving measur- 
able transformation on M^. Then 



{K{x, y) - K{x, (f){y))) uj{y) dy 



< c sup ip[x - (j){x)) (ll^ll^i + ll^ll^oo) . 

(2.3) 



Proof. Set r = sup^|x — (p{x)\; as in the proof of Lemma ||, the interesting case 
is when r < 1. We split the integral in ( |2.3|) into two parts. First, consider 



\K{x,y) ~ K{x,(j){y))\ \uj{y)\dy 

\x—y\<2r 

27r J\x-y\<2r F " Vl J \x~y\<2r F " 012/) I 

-^ (i T+i IT^ )^^II^IIl°° 

^^ J\x-y\<2r\\x-y\ F-0(y)|/ 
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If 1 

<-/ r d|/ ||u;||^oo = cr ||cj||^^ . (2.4) 

^ J\x-y\<2r F ~ y\ 

The last inequality above holds because is area preserving, and among all 
such transformations, the symmetric map = e maximizes the integral over 

Next, we consider the case when |x — y| > 2r. Observe that 

\x - </)(y)| >\x-y\-\y- 0(y)| >\x-y\-r>^\x-y\, (2.5) 

so that 

/ \K{x,y)~ K{x,(l){y))\\uj{y)\dy 

J\x-y\>2r 

'^^ J\x-y\>2r \x-y\ \x - (j){y}\ 

If r 

< - / 1 ^k(2/)|dy 

""" ''|a;-S/|>2r F ^ y\ 

^ V~'2r<|x-|/|<2 F ~ y\ J\^-y\>'^ 1"^ ~ y\ 

^ r f f^dp .. .. 1 „ „ 

<C¥.(r)(||^||^, + 11^11^^). (2.6) 

By combining the two estimates we complete the proof. D 

Finally, we give the corresponding result for the vortex method kernel. 
Lemma 5. There exists a constant C2 which is independent of a, such that 
sup\K'^{x,y) ~ K^{x',y)\ <^v(^^) . (2.7) 



Proof. Note that on R\ K'^^x, y) = K'^dx - y\) = V^G"(|a; -y\), where 

and Kq denotes the zero order modified Bessel function of the second kind ^. 
For simplicity, we take a = 1 and compute 

clC"^ 1 /l \ 1 

V) = — - - Ki{r) = — r Inr + 0{r) , (2.9) 



dr 27r V r J Air 

d^G", . 1/1 ,_ . 1 ,_ A 1 



dr^ 27r \ r^ r J An 



— - + Ko{r) + - K,{r) = - Inr + 0(1) , (2.10) 



10 
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as r -^ 0. Set r = \x — x'\ and assume, without loss of generality as K°' is 
bounded, that r < 1. 

If |a; — t/| < 2r, then \x' — y\ < \x' — x\ + \x — y\ < 3r, so that 

\K''{x,y)-K''{x',y)\ < I V^G"(|x - y|)| + |V^G"(|x - y 



< 



dG' 



dr 



i\x-y\) 



dC 



dr 



i\x-y\) 



< — r |lnr| + 0{r) . 

IT 



(2.11) 



Since dG"/dr is continuous and decays at infinity, this implies a bound of the 
form 



\K"{x,y) - K^{x' ,y)\ < cifilx - x'\) . 



(2.12) 



If, on the other hand, \x — y\ > 2r, we use the mean value theorem to 
estimate 

|K"(x,y)-K"(x',y)| < sup \VK'^{x",y)\\x-x\ 

x"eB{x,r) 

d^G" 

— {\x"-y\) r 



< 



dr^ 



< -— r Inr + 0{r) , 

4:71 



(2.13) 



which again implies a bound of the form ( p.l2| ) . In the last step we have used 
(|2.1CI|) in conjunction with \x" — y\ > r. 

To recover the scaling of the estimate in a, divide (p.l2|) by a, rescale x, x', 
and y by a~^, and note that K"{r) = K°'^^{r/a)/a. D 



C2 {X ~ X 



Corollary 6. For q G M{R^), 

f \K^{x,y) - K^{x',y)\\q{y)\dy <^-^v('- 

3. Well-posedness 



I^IIai 



:2.14) 



We can now prove the existence of unique, global, weak solutions to the 



Lagrangian flow equation (|1.9|). 

Proof of Theorem [J. Due to the quasi-Lipschitz condition for if", we can adopt 
the method that Kato developed for the Euler equations in [|I^, by simply re- 
placing the kernel estimates in L^ by the corresponding estimates in L°°. Our 



presentation follows to some extent that of Marchioro and Pulvirenti [ pO[] . 
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For simplicity, we assume a = 1 throughout this proof. We introduce a 
sequence of approximate solutions 

dtr]''{x,t)=u''{r]''{x,t),t), (3.1a) 

ri''{x,0)=x, (3.1b) 

r]\x,t) = x, (3.1c) 

q^{r{x,t),t)=qo{x), (3.1d) 

«"(x, t)= [ K"(x, y) q--\y, t) dy , (3.1e) 

for n eN. The proof now proceeds in several steps. 

Step 1. Prove that ?]" G C^((0, oo); Q) for every n eN. 

We proceed inductively. Notice that for every n the vector field u"' is quasi- 
Lipschitz in space and continuous in time. This is a consequence of Lemma ^ 
as 



|M"(x,t)-M"(a;',t)| 



K"{x, v"-\y, t)) - K%x\ v^-\y, t)) qo{y) dy 



< Cip{x — x') \\qo\ 



M 1 



(3.2) 



and 



k"(x,t)-M"(a;,t') 



K^{x,i^'^-\y,t)) - K^{x,7^'^-\y,t'))\q,{y)dy 

< sn^\K"{x,^^-\y,t)) - K''{x,t-\y,t'))\\\qo\\M 
<c sup if {r]''-\y,t)-r]^-\y,t')) ||golU 



< c sup sup v?(|^" ^iy^s)\\t-t'\)\\qo\\j^ 

yeR^ x(^[t,t'] 

= c sup sup v5(|m""^(2/, s)| |t-t'|) llgollx 

" xG[t,t'] 



(3.3) 



This implies uniform continuity in time, because m" is bounded for every n: 



K{x,t)\ 



K-ix,r]^-\y,t))qoiy)dy 



< snp\K{x,y)\ ||go|U = c|kolU • 

(3.4) 

Since m" is continuous in time and quasi-Lipschitz in space, the vector field 
generates a local flow i]"' G C^([0,T); C(M^)) for some T > — see, e.g.. Chap- 
ter 2, Lemma 3.2 in Marchioro and Pulvirenti |^. Because of the global 
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bound (p^), the right side of (|3.1a|) is bounded and the flow exists globally in 
time. 

Step 2. Show that there exists a limiting flow map rj G C((0, oo); Q). 

We flrst prove that the sequence i]^ is Cauchy in C([0, T]; Q) for some T > 0. 
To simplify notation, we shall drop the explicit time dependence of u and rj, 
and estimate 

\r{x,t)-r'\x,t)\ 

< [ |M"(r7")-u"-^(r7"-^)|ds 



< 



ir"(r/"(x),r^"-^(y))-7r"(r^"-i(a;),r/"^^(y)) qoiy) dy 



ds 



+ 



ir"(r/"-i(x),r/"-^(l/))-ir-(r/"-i(x),r/"-^(l/)) qo{y) dy 



ds 



<c / ip[r]"-{x) - if -^(x)) ds llgol 



A^ 



+ c y.(r/"-i(x)-r^"-2(x))ds||golU- 
By taking the supremum over x on both sides, we obtain 

s\x^\rf{x,t) - 7f-^{x,t)\ 



(3.5) 



< c||go| 



M 



^ 



I sup|?7"(x)-?7" ^(s)| ) +(/?( suplr/" ^(s)-//"" 



Defining 



p^^{t) = sup sup|?7"(x, t) -7]'"' {x,t)\ 

n>N xGR'^ 

we can simplify the previous estimate, and obtain 

p^(t)<c/V(p^-i(s))ds. 



It is well known that this implies 



limp^^(t)^0 

AT— >oo 



ds. 
(3.6) 

(3.7) 

(3.8) 
(3.9) 



uniformly on [0, T] for T sufficiently small. Since T depends only on a and the 
7W-norm of qo, this result can be extended to arbitrarily large times. Thus, 
the contraction mapping theorem implies the assertion of Step ^ 
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Step 3. Show that the Lagrangian flow equation ( |1.9|) is satisfied in the hmit, 
and that r] e C\R;g). 

We define the hmiting potential vorticity q and the limiting velocity u in 
the obvious way, and check by direct estimation that 

q"^q = qoor]-' (3.10) 

weakly in A^(M^), and 

u''^u = K''*q (3.11) 

in C(R^); both limits are uniform over finite intervals of time. 

To prove that i], u, and q solve the limit problem (p..9|), we consider its 
integrated version 

r]{x^t) — / u{ri{x, s), s)ds 
Jo 

= ■r]{x,t) — / u{ri{x,s), s) ds — rj"'{x,t) + / u^{ri^{x,s),s) ds 
Jo Jo 

r\ 

< |?7(x, t) - ?7"(x, t)| + / |M"(r7"(x, s),s) - ■u(?7"(x, s),s)| ds 

Jo 

[\ I 

+ / \u{rf'{x,s),s) — u{r]{x,s),s)\ds 
Jo 

—>■ uniformly in a; as n — > oo . (3.12) 

Thus, the left side must be zero. Since u{ri{x,s),s) is continuous in x, we 
can differentiate with respect to t, and find that t] satisfies ( |1.9| ) and that i) 
is in fact continuous. Due to the time-reversibility of the equation, the result 
extends to negative times as well. 

Moreover, one can show — first by formal calculation for smooth function, 
and then extending by the usual density argument — that the weak solution q 
defined through ( ^.10| ) satisfies 

/ / ((9i0 + M-grad0)gda;dt = (3.13) 

Jr Jr2 

for every G C^(]R x R^). This shows that solutions of the vortex method, 
and hence the equations of second-grade non-Newtonian fiuids, preserve the 
(weak) co-adjoint action. 

Step 4. Prove that the solution is unique. 
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Uniqueness is shown by a direct estimate on the difference of two fiow maps. 
This leads to another log-Gronwall inequahty, which can be treated in the 
same way as the previous ones; we omit all details. D 



Remark 8. The homeomorphisms that we consider have the vector space M? 
as the range; we may thus subtract two elements of this class. For homeo- 
morphisms of a compact domain VL of M^ , one can isometrically embed the set 
of measure-preserving homeomorphisms of Vt into the vector space L^(r2,R^), 
and take differences in this large space. Similarly, the difference m° o r/" — m o r/ 
is not an intrinsic operation, but rather relies on the trivial identification of 
vector spaces induced by the trivial geometry of M^. On the other hand, when 
the configuration space is ©^(M), s > 2 and M is a compact Riemannian 
manifold, the map u"" o r^" is an element of the fiber T^aV'l while m o r^ is in 
TrfDl^\ thus, in order to compare the two maps, we must parallel transport uor] 
into TrjcVf^ along the Riemannian connection. 

4. Weak co-adjoint action and reduction 

As we described, classical solutions of the two-dimensional averaged Euler 
equations are geodesies on the Hilbert-class volume-preserving diffeomorphism 
group P^, s > 2. We identify the space of classical vorticity solutions H^~^{M) 
with the reduced space TeP^ = TP^/P^ (symmetry reduction by the massive 
particle relabeling symmetry group P^ of hydrodynamics), and note that this 
space is the union of the P^-co-adjoint orbits. 

In the case that M = M?, and for the purpose of studying weak solutions to 



(|1.2| ) we shall substantially relax the regularity requirements on the configura- 
tion space, and use Q in place of P^; correspondingly, we shall use the vector 
space of Radon measure on R^, which we denote by A^(R^), for the reduced 
space of vorticity functions, in place of the space of if ^~^ functions. 
Recall that the co-Adjoint action of P^ on if^~^(R^) is given by 



M 



Ad;(g) = go77. (4.1) 

We shall need to define the notion of weak co- adjoint action of Q on M.{E?). 
First, note that the operation Ad* : Q x M{^^) ^ M{^^) given by Ad*(g) = 
q o Tj is well-defined. Next, define the weak co- Adjoint action w-Ad* : Q x 
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7V1(M2) ^A<(r2) by 

/ / w-Ad*^{q)-(j)dxdt= I I q ■ {(j) o 7]) dx dt (4.2) 

for all G C^(R x M^). 

It follows that if rjt is a C^ curve in ^ such that e = rjQ and m = {d/dt)\t=oi]t, 
then we may — computing the time derivative of w-Ad* (g) at t = — define 
the weak analogue of the algebra co-adjoint action by 

/ / w-ad*(g) ■ 0dxdt = / / g ■ ((9t0 + m ■ grad (/)) da; dt (4.3) 

Jr J r2 Jr iiR2 

for every G C^(]R x M^). Recall that the classical co-adjoint action is defined 
by ad*Jgi) = {d/dt)\t=o{r]lqt) where {d/dt)\t=o{r]tQt) = dtqt + C^St- In two 
dimensions, the Lie derivative term CuSt reduces to Ut ■ gradg*. 

Theorem 7. For any go ^ A^(M^), let Oq^ denote the co- Adjoint orbit {q'- q = 
qooi],!] G Q}. The weak co-adjoint action o/C^j^(M^) on A^(M^) is well-defined, 
and solutions of the second-grade fluids equations or of Chorin's vortex blob 
method with initial data go leave Oq^ invariant. 

Proof. The result immediately follows from the fact that the vanishing of the 
weak co-adjoint action is equivalent to the weak formulation of ( |1.2| ). The- 
orem [^, giving global well-posedness of weak solutions, then concludes the 
argument. D 

5. Convergence 

We can now prove convergence of the flow of the vortex blob method to the 
flow of the Euler equations. This is done in two steps. First we show that 
the averaged Euler equation, or vortex method PDE, approximates the Euler 
equation as a — > for bounded vorticity fields. In the second step, we prove 
that continuous solutions of the averaged Euler equation can be approximated 
by measure- valued ones. These two results together imply Theorem |^. 

Lemma 8. Let go = cuo G L^(R^)nL°°(R^). Then for every T > there exists 
a positive constant C{T) such that 

sup snp Iri" {x, t) - ri{x, t)\ < C (T) a""'^ . (5.1) 

te[0,T]iGK2 
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Proof. We estimate the difference of the Euler and Euler-a ffow maps: 
\ri"{x,t) — ri{x,t)\ < / \u" o T]" — u o rjl ds 



JO 

< [ [ \K'^ir^'^ix),v^iy)) - Kir^'^ix),v''{y))\\coo{y)\dyds 

Jo JK2 

+ / / \K{r]"{x),r]"{y))-K{r]ix),r]^{y))\\uJoiy)\dyds 
Jo Jm.'^ 

+ \K{r]{x),r]''{y))-K{r]{x),ri{y))\\uJo{y)\dyds 

Jo Jr'2 
ft 

{I, + I^ + I,)ds (5.2) 

'0 

To estimate /i, we note that on M^, the difference of the kernels is exphcitly 
given by K°'{r) — K{r) = Ki{r/a)/a, so that 



h 



-Kj\^^^]\u;^iy,s)\dy 
]R2 a \ a 



< ca llcjoll^oo (5.3) 

The other two integrals can be estimated by using the quasi-Lipschitz condi- 
tions, Lemma |^ and Lemma ^, respectively. One finds that 

h < Cip{r]°'{x, s) - r]{x, s)) (Hc^oll^i + IkoH^oo) , (5.4) 

and 



h= \K{ri{x,s),y)- K{ri{x,s),ri°' OT] ^{y, s))\\uj{y, s)\dy 

JR2 

< c snpip{x-r]''or]'\x,s)) (||^(s)||^i + \\uj{s)\\^^) 
= c sup ip{ri{x, s) - ?7°(x, s)) {\\ujo\\^i + H^oll^oo) . 



(5.5) 



By inserting the bounds for Ji to Is back into (|5.2| ) and taking the supremum 
on both sides, we obtain the log-Gronwall inequality 

sup |r7°(x, t) — ?7(x, t)| < / a Ki + K2 snp (p(^ri°'{x, s) — ri{x, s)) ds . (5.6) 

XgM2 Jq I ^g]g2 J 

To obtain explicit bounds that are valid on any finite interval of time [0,T], 
we set 



p{t) = sup \ri°'{x, t) — ■r]{x, t) I 



(5.7) 
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and use the tangent approximation of the concave function 99; namely, for any 

£€(0,1), 

ip{r) < ip{e) + ip'{e) r = {-\ne)r + e . (5.8) 

This makes the right-hand-side of ( |5.8| ) hnear in r. For notational simphcity, 
we also rescale a and t such that Ki = K2 = I. We substitute ( |5.8| ) into (^]B|) 
and obtain the usual Gronwall inequality; it follows that p must satisfy the 
differential inequality 

p<{-lne)p + e + a, p(0) = . (5.9) 

Setting e = e'^a"^^^^'^^ and integrating (|5.9| ) with this choice of e{a), we find 
that 

pit) < a^ + e* . 5. 10 

e — ma 

Thus, p = ©(a'^^P^-^)) uniformly on [0, T]. D 

In the following we will consider a as fixed and approximate continuous data 
by measure valued data. Let 77, q, and u denote quantities corresponding to 
a solution of the Euler-a equation with initial data go ^ L°°(R^), and let 77", 
g", and m" denote a sequence of solutions to the Euler-a equation with initial 
data Qq G 7W(]R^) for every n G N. Then the following is true. 

Lemma 9. Let qo G L (R^) fl L°°(R^), and suppose that go ^-5 approximated 
by a sequence of measures in 7V1(R^) such that q^ -^ go weakly in M.{E?), and 
WIqWm ~^ Ikollii- Then, for every T > 0, 

lim sup sup |?7"(x, t) — 77(x, t)| = . (5-11) 

"^°°ie[o,T]xeiR2 

Proof. As in the proof of Lemma p, we estimate 
|r/(x,t)-r/"(x,t)| 



<-f 

Jo 

<-f 

Jo 



[K'^ivix, s),viy, s)) qoiy) - ir"(r7"(x, s), r^"^/, s)) q^iy)] dy 



ds 



K"{7l{x),r]{y)){qo{y)-q^{y))dy 



ds 



+ [ [ \K'^{r]ix),r]iy)) - K'^{r]ix),v"iy))\\qoiy)\dyds 

Jo JR2 

+ r / \K^{v{x),v''{y)) - K^{v"{x),v''{y))\\qo{y)\dyds 

Jo ^R2 
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= / iJl + J2 + J3)ds 

Jo 
We find, after a cliange of variables, tliat 



sup Ji = sup 



K^{x,y){qo-q^){v"\y))dy 



(5.12) 



(5.13) 



By Lemma p!0| below with 0(x ~ y) = K°'{x,y) and qn{y) = {qo — <lo)iv ^iv)): 
this expression converges to zero as n ^ oo. Moreover, by Lemma |^, 

J2< sup snp\K'^{rj{x),rj{y))-K'^{r^{x),r]^{y))\U\\M 



<sup^^(^M^^^')|KI 



and 



a 



c f r]{x)-rj''{x) \ 
h<-^{ ]ho\\M 



M 



(5.14) 



(5.15) 



By inserting these estimates back into (|5.12| ) and taking the supremum in x 
on both sides, we obtain an integral inequality that can be solved with the 
log-Gronwall inequality exactly as in the proof of Lemma ^ The result then 
follows. n 

Lemma 10. Let {g„} be a sequence of measures in A^(R^) converging weakly 
to zero with uniformly bounded total variation and uniform decay at infinity. 
Further assume that (p is a continuous test function with — i> as |a;| — > oo. 
Then 



lim sup 



{x - y) qn{y) dy = 



(5.16) 



Proof. Set M = sup„||grt||A^ and M' = sup^\(f){x)\. Let e > be fixed. By 
assumption on the {g„}, there exists an _R > such that for every n G N, 

e 



\(lniy)\dy< 



\y\>R 



6M' 



(5.17) 



Moreover, there exists a.n R' > such that \(p{x)\ < e/(2M) for |a;| > R'. 
Since is uniformly continuous on compact sets, there exists 6 > such that 
|0(x) -0(x')| < e/{3M) for allx,x' G B{0,2R + R') with \x - x'\ < S. Cover 
B{0,R + R') with finitely many balls of radius 6 and denote the centers of 
these balls by Xj, i G /. Choose A^ large enough such that for n > N, 

e 



max 

i&I 



[Xi - y) qniy) dy 



< 



(5.18) 
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Then for |a;| < i? + R' there exists an i G / such that |a; — XjI < 5, and 
[x - y) qniy) dy 

(j){x-y)-(f){xi-y)\\qn{y)\dy+ / (j){xi - y) qniy) dy 

< sup|0(x-|/) -0(a;i-|/)| ||g„||_^ + 2sup|0(a;)| / \qn{y)\dy . 
\y\<R xm? J\y\>R "^ 

< ^M + 2M'^ + - = e. (5.19) 
- 3M QM' 3 ^ ^ 

On the other hand, if |x| > R + R', then 

Kx-y)qniy)dy 



+ 



< sup 10 (a:) I / \qn{y)\dy+ snp\(f){x)\ / \qn{y)\dy 
\x\>R' J\y\<R ^eM2 J\y\<R 



< —M + M' ^ 



2M 
This completes the proof. 



6M' 



< e . 



(5.20) 
D 
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